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I. 


INTRODUCTION 


The  long  range  goal  of  this  research  is  to  develop  techniques 
for  detecting  and  identifying  (or  classifying)  ship  targets  using  H.F. 

Over  the  Horizon  radar  (OTH).  For  the  period  covered  in  this  report, 

(the  initial  phase)  several  specific  studies  were  undertaken.  The 
general  identif iabi lity  of  naval  surface  ships  using  H.F.  radar  was 
studied  first.  This  was  a  study  of  the  behavior  of  an  available  data 
set  from  the  point  of  view  of  radar  target  classification.  In  this 
sense,  it  was  not  restricted  to  a  study  of  the  performance  of  any  specific 
target  classification  algorithm.  It  was,  rather,  a  quantitative  appraisal 
of  the  overall  usefulness  of  H.F.  radar  returns  for  target  classification. 
Two  techniques  were  used  for  this  identif iabi 1 ity  study.  The  nearest 
neighbor  technique  was  used  to  study  the  clustering  property  of  the 
data  in  multi-dimensional  space,  while  the  concept  of  the  displacement 
and  skew  of  lines  in  multidimensional  space  (hyperlines)  was  used  to 
study  the  separability  of  the  data.  Time  domain  and  resonant  pole 
techniques  were  also  implemented  and  applied  to  the  available  data 
base. 


The  results  of  these  studies  indicate  that  data  from  H.F.  radar 
systems  are  useful  for  ship  target  classification.  In  the  language 
of  pattern  recognition,  the  data  were  shown  to  cluster  in  hyperspace 
and  the  clusters  were  shown  to  be  separable.  The  time  domain  studies 
were  used  to  show  that  the  main  contributors  to  the  backscatter  features 
are  large  significant  ship  structures  such  as  the  bridge  geometry. 

A  quantitative  measure  of  these  statements  will  be  given. 

II.  DATA  BASE 

The  automatic  identification  or  classification  of  an  unknown 
member  (ship)  of  a  class  (of  ships)  based  on  various  subsets  of  a  set. 
of  features  (measured  cross  section  and  phase  vs.  H.F.  frequency)  was 
studied.  The  first  study  undertaken  was  to  measure  the  inherent  identi- 
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fiability  of  ships  based  on  H.F.  OTH  radar  returns.  Once  this  property 
has  been  answered  quantitatively,  the  development  of  specific  applicable 
algorithms  and  their  optimization  can  proceed.  Both  classifiability 
tests  and  the  testing  of  specified  classification  algorithms  require 
a  data  base  which  is  representative  of  the  data  to  be  expected  in 
practice.  This  classif iabi lity  study  can  work  with  a  relatively  small 
data  set. 

A  set  of  radar  backscatter  measurements  made  by  the  Ohio  State 
University  ElectroScience  Laboratory  of  scale  models  of  ships  was  used 
as  the  data  base  for  the  studies  reported  here  [1,2,3].  A  summary  of 
the  data  base  is  shown  in  Table  1.  In  each  case,  scale  models  of  the 


TABLE  1 

OSU  DATA  BASE  (1979) 


SHIP 

SCALE 

TYPE 

SCALED  Fn 
(mHz)  U 

#  OF  HARMONICS 

Midway 

(1 : 500) 

Carrier 

2.17 

10 

Missouri 

(1:500) 

Battleship 

2.17 

10 

Missouri 

(1:700) 

Battleship 

1.55 

10 

Sverdlov 

(1:500) 

Cruiser 

2.17 

10 

B ismark 

(1 :700) 

Battl eship 

1.55 

10 

Mogami 

(1:700) 

Heavy  Cruiser 

1.55 

10 

Hayanami 

(1:700) 

Destroyer 

1.55 

10 

Shimokaze  (1 :700) 

Destroyer 

1.55 

10 

All  ships  measured  at  bow,  stern,  port,  and  starboard  aspect 
angles. 

The  Missouri  (1:700)  and  the  Sverdlov  were  measured  at  5° 
aspect  angle  increments  about  the  bow,  stern,  and  starboard 
aspects. 


listed  ships  were  measured  with  a  29°  elevation  angle  using  a  multi  frequency 
bistatic  radar  system  (1).  The  models  were  placed  on  a  conducting 
ground  plane.  Measurements  at  5°  aspect  angle  increments  were  made 
for  the  Sverdlov  and  the  Missouri  as  shown.  The  scaled  fundamental 
frequency  is  shown  in  the  table.  Ten  harmonically  related  measurements 
were  made  of  the  absolute  cross  section  "and  the  absolute  phase  (relative— 
to  a  point  defined  at  the  center  of  the  hull  of  the  ship).  As  can 
be  seen,  these  measurements  (scaled)  cover  a  large  part  of  the  H.F. 
band.  A  photograph  of  these  ship  models  is  given  (Photograph  1). 

It  is  important  to  note  that  these  ships  are  all  generally  of  the  same 
size,  and  that  the  size  of  a  particular  ship  will  not  be  the  determining 
characteristic  in  the  tests  to  be  performed  on  the  data. 


The  data  for  the  Missouri  (1:700)  and  the  Sverdlov  were  used 
for  the  identif iabi 1 ity  testing  since  a  much  larger  set  of  aspect  angles 
was  available  for  them.  An  example  set  of  data  comprising  one  frequency 
.  sr.au  i s  sh^wo  .i n..F j gure_J . . .  Jhi s  ^phasor,  rep r  ese nt a ti o  n of  the 

measured  amplitude  square  root  of  (fcross  section '( cm) )  and  absolute  phase 
of  the  Sverdlov  at  an  aspect  angle  of  5°.  The  entire  data  base  is 
given  in  Appendix  A  in  this  compact  phasor  plot  form.  Study  of  the  data 
in  the  appendix  have  shown  several  interesting  properties. 


First  note  that  as  the  aspect  angle  of  a  particular  ship  changes  by 
5  nr  10  degrees,  the  position  in  phasor  space  of  harmonic  N  (N  =  1  to  5) 
remains  relatively  fixed  (compared  to  the  corresponding  phasor  position  of 
the  data  from  the  second  ship).  For  both  ships,  the  broadside  amplitudes 
tend  to  be  3:  4  times  larger  than  the  bow  and  stern  amplitudes.  Distinct 
differences  in  the  behavior  of  the  Sverdlov  data  as  compared  with  the 
Missouri  data  can  also  be  seen  in  the  phasor  plots.  The  techniques  used  in 
the  next  three  sections  of  this  report  can  be  used  to  quantify  this  obser¬ 
vation. 
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figure  1.  Example  measurement  (Sverdlov  at  5°  off  bow).  Phaser  plot 
of  the  amplitude  (square  root  of  the  radar  cross  section  - 
measured  in  cm?)  and  phase  of  the  data  for  the  Sverdlov. 
Numbers  are  the  harmonic  multiples  of  2.17  MHz. 


III.  NEAREST  NEIGHBOR  TECHNIQUES 


The  data  from  the  H.F.  radar  return  of  a  ship  can  be  represented 
by  a  \<>t  of  points  in  a  multidimensional  space.  Each  measurement  parameter 
is  represented  in  a  different  dimension  of  the  space.  Each  measurement 
at  a  set  of  frequencies  on  a  ship  will  yield  a  point  in  this  multi¬ 
dimensional  space  (hyperspace).  Measurements  at  a  number  of  elevation 
and  aspect  angles  will  yield  a  set  of  such  points.  The  behavior  of 
the  set  of  points  determines  the  overall  ability  of  the  radar  system  to 
identify  a  target.  It  is  necessary  for  the  data  processor  of  the  H.F. 
radar  system  to  decide  which  of  several  candidate  groups  of  points 
a  new  unknown  point  belongs  to. 

A  number  of  properties  of  a  set  of  such  points  in  hyperspace 
can  be  defined.  For  example,  the  points  for  a  particular  ship  may 
all  lie  in  a  localized  region  of  hyperspace  (i.e.,  they  may  exhibit 
clustering).  This  would  be  a  useful  property.  Alternatively,  the 
points  may  lie  on  a  line  or  a  surface  (actually  a  hyperl i ne  or  a"yhyper^* 
surface).  Conceptually,  it  is  clear  that  it  is  not  sufficient  for 
the  points  measured  from  a  particular  ship  to  cluster  in  order  for 
the  classification  system  to  be  successful.  It  is  also  necessary  for 
the  data  representing  the  other  ships  in  the  study  to  occupy  some  other 
location  in  hyperspace.  Nearest  neighbor  techniques  [4]  will  be  used 
to  study  clustering  of  the  ship  data  in  this  report,  while  a  hyperline 
technique  will  study  the  property  of  separation  of  the  data  into  distinct 
regions  of  hyperspace. 

The  technique  used  to  study  the  clustering  property  of  the  data 
base  available  here  is  called  the  nearest  neighbor  technique.  This 
technique  is  non-parametric  and  therefore  the  statistical  properties 
of  the  data  need  not  be  known  for  the  results  to  be  valid.  The  concept 
of  the  nearest  neighbor  calculation  is  simple.  Each  set  of  measure¬ 
ments  of  the  ship  target  yields  as  many  as  20  distinct  parameters 
(i.e.,  the  real  and  the  imaginary  parts  of  the  radar  cross  section 
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it  each  of  t. ho  10  frequencies).  This  means  that  each  aspect  angle 
measured  for  the  ship  yields  a  point  (or  vector)  in  multidimensional 
space  and  that  the  space  may  be  defined  by  as  many  as  20  dimensions. 

The  distance  between  a  particular  point  and  each  of  the  other  points 
in  this  space  may  be  computed,  and  the  smallest  of  the  distances  is 
taken  as  the  distance  to  the  nearest  neighbor  of  the  original  point. 

The  original  point  may  represent  data  from  an  unknown  target  and  the 
other  points  represent  data  from  a  set  of  cataloged  measurements  of 
known  targets  at  various  aspect  angles  and  elevation  angles.  The  unknown 
point  is  associated  with  the  target  from  which  the  nearest  neighbor 
was  derived.  This  concept  forms  the  basis  for  a  target  classif ication 
algorithm. 

In  this  report,  the  nearest  neighbor  technique  will  be  used  to 
test  the  clustering  property  of  the  data  base  rather  than  strictly 
as  a  target  identification  algorithm.  In  this  way,  a  more  fundamental 
property  of  the  original  data  set  may  be  assessed.  The  data  from  a 
single  frequency  scan  of  a  target  (the  Missouri  (1:700)  or  the  Sverdlov 
-  as  mentioned)  for  a  subset  of  the  harmonics  available  and  both  with 
and  without  phase  data  are  used  to  derive  a  set  of  points  (vectors) 
in  hyperspace.  Then  the  distance  from  the  point  in  question  to  all 
of  the  other  points  was  computed.  The  minimum  of  these  distances  was 
used  to  define  the  nearest  neighbor.  For  this  report,  a  set  of  ship 
data  is  clustered  if  the  nearest  neighbor  to  a  particular  point  (or 
vector)  is  a  point  from  the  same  ship  (at  some  other  aspect  angle). 

An  example  of  the  results  from  this  study  is  given  in  Table  2. 
Abbreviations  are  used  to  identify  a  particular  set  of  data.  The 
Missouri  data  measured  at  an  aspect  angle  of  10°,  for  example,  is  listed 
as  "MISOIO".  The  left  hand  column  lists  the  test  data  (or  unknown 
data),  and  the  candidate  data  used  for  comparison  testing  is  listed 
at  the  top  of  each  column.  The  numbers  are  the  actual  distances  which 
were  computed;  they  have  not  been  normalized.  At  the  top  of  the  first 
column  is  the  first  harmonic  tested  (1)  and  the  last  harmonic  tested 


TABLE  2 


(10).  The  T  indicates  that  phase  measurements  were  used  to  convert 
the  data  to  real  and  imaginary  components  prior  to  computing  the  neighbor 
distances.  There  were  a  total  of  20  dimensions  used  in  the  example 
shown.  The  example  shown  has  had  the  nearest  neighbor  value  (minimum) 
underlined.  If  the  nearest  neighbor  is  in  fact,  the  other  ship,  the 
value  is  underlined  twice.  The  entire  set  of  nearest  neighbor  results 
is  given  in  Appendix  B.  The  data  given  in  Appendix  B  (and  in  Table  2) 
give  a_M  distances  in  multidimensional  spate  between  the  test  target 
and  the  catalogued  data  -  no  normalization  has  been  used.  A  total 
view  of  the  clustering  properties  of  the  data  base  can  be  found  in 
this  appendix.  Some  of  the  general  behavior  of  the  nearest  neighbor 
calculations  for  all  the  data  can  be  seen  in  the  example  shown  in  Table 
?.  It  can  be  seen  that  the  nearest  neighbor  is  typically  not  only 
from  the  same  ship  as  the  test  (unknown)  data,  but  is  also  close  in 
aspect  angle  to  the  test  data. 

Note  that  the  comparisons  between  identical  data  sets  will  always 
produce  a  zero  vafue  in  this  noiseless  or  errorless  type  of  study. 

If  the  statistics  of  the  errors  can  be  defined  in  a  meaningful  way, 
then  artificially  introduced  signal  corruption  could  be  included  in 
the  study  and  a  measure  of  the  expected  performance  of  the  nearest 
neighbor  algorithm  for  target  identification  could  be  studied.  This 
type  of  study  will  be  included  in  future  research  tasks. 

The  example  shown  here  has  an  anomalous  point.  The  nearest  neighbor 
of  the  Sverdlov  at  175°  aspect  angle  (SVE175)  is  a  data  point  for  the 
Missouri  at  5°  aspect  angle  (MIS005).  Note  however,  that  the  nearest 
neighbor  in  this  case  is  not  very  near.  In  fact,  the  nearest  neighbor 
distance  for  SVE175  is  more  than  an  order  of  magnitude  larger  than 
several  of  the  other  nearest  neighbor  distances.  Note  (assuming  the 
data  base  is  accurate  here)  that  this  error  does  not  mean  that  the 
nearest  neighbor  test  for  target  identification  would  necessarily  fail 
in  this  case.  It  is  likely  in  practice  that  the  catalog  of  data  would 
contain  measurements  of  the  Sverdlov  at  175°  for  comparison  with  an 
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unknown  target  which  would  also  be  the  Sverdlov  at  175°  in  this  case.  Only 
the  errors  in  the  measurements  would  contribute  to  miscl assif ication. 

On  the  other  hand,  if  the  data  base  is  valid,  this  test  has  identified 
a  sensitive  area  for  further  study.  It  can  be  seen  that  small  changes 
in  the  aspect  angle  have  introduced  large  changes  in  the  features 
of  the  H.F.  radar  return,  namely  the  data  are  not  clustered  at  this 
aspect  angle. 


A  summary  of  the  nearest  neighbor  tests  is  given  in  Table  3. 


TABLE  3 
NN  SUMMARY 


NEAREST  NEIGHBOR 


HARMONICS 

BROADSIDE 

BOW-STERN 

ALL 

ASPECTS 

AMP 

ALL  ASPECTS 

AMP 

AMP 

AMP+Ph 

+/-50 

AMP 

AMP+Ph 

+/- 5° 

P:  Q 

P:Q  ___ 

P :  Q 

P:Q 

P:Q 

P:Q 

1-5 

9:10 

9:10 

139 

139 

12:12 

3-7 

9:10 

10:10 

mm 

12:12 

6-10 

6:10 

8:10 

5:10 

12:12 

11:12 

12:12 

1-10 

5:10 

7:10 

6:10 

12:12 

11:12 

12:12 

P  =  Number  of  times  that  the  nearest  neighbor  of  a  target 
was  the  same  target  but  at  a  different  aspect  angle. 

0  =  Number  of  cases. 

P  r  Q  is  a  measure  of  the  "clustering"  property  of  the  data. 
Ties  counted  as  clustered, 

CONCLUSIONS: 

The  BOW-STERN  data  are  all  well  clustered.  Identification 
should  be  possible  even  with  large  measurement  errors. 


The  lower  order  harmonics  exhibit  more  clustering  than  the 
higher  order  for  the  broadside  data. 


Iri  this  table,  the  data  were  grouped  in  various  subsets  of  the  total 
number  of  available  harmonics.  The  ratio  of  clustered  points  (P)  to 
l.li"  number  of  comparisons  (0)  is  given  in  this  table.  The  first  three 
I'ntries  in  the  table  all  use  !>  harmonics.  This  represents  b  dimensions 
for  the  amplitude  only  cases  and  10  dimensions  for  the  amplitude  and 
phase  (real  and  imaginary  components)  cases.  Note  that  all  of  the 
bow-stern  amplitude  only  cases  are  well  clustered  (i.e.,  12:12).  A 
test  of  data  comparisons  using  only  +/-50  aspect  angles  is  also  presented 
in  this  table.  This  is  a  test  of  the  clustering  of  the  data  if  a  priori 
information  on  the  aspect  angle  is  available.  It  is  possible  to  see 
that  poor  clustering  is  not  caused  by  including  too  many  aspect  angles 
in  the  test  data.  It  can  also  be  seen  that  the  clustering  of  the  broad¬ 
side  data  is  not  improved  by  the  introduction  of  phase  data  into  the 
algorithm.  In  fact,  in  some  cases,  the  clustering  is  degraded  by  the 
introduction  of  phase  information.  Finally,  note  that  the  data  from 
the  lower  order  harmonics  exhibit  more  clustering  than  the  higher  order 
harmonic  data. 

It  should  be  emphasized  that  degradation  of  the  clustering  property 
of  a  set  of  data  does  not  mean  that  the  identification  potential  of 
the  data  has  been  reduced.  It  may,  in  fact,  have  been  increased  since 
loss  clustering  means  the  data  are  more  spread  out  in  feature  (mult:- 
dimcns ional )  space  and  thus  less  sensitive  to  noise,  clutter,  etc. 

On  the  other  hand,  reduced  clustering  as  defined  in  this  report  does 
moan  that  the  aspect  angle  sensitivity  has  been  increased  (as  might 
be  expected  for  higher  harmonic  numbers  or  the  introduction  of  phase!. 
Thus  the  catalog  of  measurements  required  for  identification  of  a 
particular  target  must  include  a  larger  set  of  aspect  (and  elevation' 
angles  if  the  clustering  is  reduced. 

In  conclusion,  then,  the  data  are  well  clustered  for  the  bow- 
■» fern  group,  and  for  the  lower  order  harmonic  values  at  broadside. 

A  priori  information  does  not  increase  the  clustering,  and  neither 
does  the  introduction  of  phase  information.  These  general  results 
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are  expected  to  be  applicable  to  the  full  scale  OTH  measurements  taken 
by  the  NRL  SEA  ECHO  system.  In  other  words,  it  is  expected  that  ampli¬ 
tude  data  alone  on  full  scale  ships  will  cluster  well  for  bow-stern 
and  will  probably  cluster  somewhat  at  broadside  aspect. 

IV.  HYPERLINE  STUDIES 

It  is  not  sufficient  for  data  in  multidimensional  space  to  cluster 
in  order  to  be  able  to  confidently  predict  the  identif iabi 1 ity  of  the 
ship  targets  using  H.F.  radar  return  information.  The  clusters  must 
also  separate  from  each  other  by  an  amount  greater  than  the  expected 
measurement  errors.  Algorithms  for  exploiting  the  data  for  target 
identification  may  then  be  confidently  developed  and  tested.  Furthermore, 
clustering  is  not  the  only  property  that  predicts  the  identif i abi 1 ity 
of  the  ship  targets  in  the  data  base.  It  can  be  shown,  for  example, 
that  the  data  points  in  multidimensional  space  derived  in  the  study 
in  this  report  should  form  closed  curves.  This  may  be  seen  if  it  is 
noted  that  the  only  variable  in  the  data  base  presently  in  use  is  the 
aspect  angle.  The  only  geometric  configuration  that  results  from  a 
function  with  a  single  independent  variable  (modulo  ?ir)  is  a  closed 
curve.  This  curve  may  be  quite  complex,  of  course,  but  this  complexity 
is  constrained  by  the  observation  from  the  data  base  used  here  that 
there  is  a  reasonably  smooth  progression  in  the  backscatter  data  as 
the  aspect  angle  is  changed.  There  are  also  theoretical  limits  to 
the  complexity  of  these  closed  curves.  These  are  related  to  the  size 
of  the  ships  in  proportion  to  the  wavelength  at  H.F.  (See  the  section 
on  time  domain  studies.) 

The  study  undertaken  to  quantify  these  ideas  was  taken  from  a 
representation  of  the  data  as  points  (or  vectors)  in  a  multidimensional 
space.  The  curves  which  result  from  a  variation  in  aspect  angle  for 
the  two  targets  will  be  separated  from  each  other  by  some  distance, 

D.  this  distance  may  be  computed  for  the  given  data  over  restricted 
ranges  of  aspect  angles.  The  method  is  illustrated  by  Figure  2.  It 
is  summarized  below. 


A  DATA  SET  FOR  TARGET  I 
O  DATA  SET  FOR  TARGET  2 


Figure  2.  Hyperline  concept  -  definitions  of  D  and  *. 


(1)  Restrict  the  range  of  aspect  angles  (and  number  of  harmonics). 

(2)  Fit  a  straight  line  to  each  set  of  points  over  the  range  of  aspect 
angles  chosen.  (The  "goodness  of  fit"  of  the  set  of  points  to 

a  straight  line  will  determine  the  applicability  of  the  data 
to  this  type  of  treatment.) 

(3)  Compute  the  distance  D  between  the  two  straight  lines  in  multi¬ 
dimensional  space  (hyperlines)  at  their  closest  approach,  and 
the  skew  angle,  G,  between  the  two  hyperlines.  The  computations 
used  and  the  derivation  of  the  equations  are  given  in  Appendix 
C. 

This  procedure  reduces  the  description  of  the  separation  of  two 
mult i parameter  feature  sets  to  a  scalar  number  D.  The  value  of  D  can 
be  computed  as  a  function  of  the  number  of  frequencies  used  and  the 
range  of  aspect  angles  chosen.  The  parameter  0  is  a  measure  of  the 
skewness  of  the  two  hyperlines.  This  parameter  gives  a  measure  of 
the  sensitivity  of  the  distance  D  to  changes  in  aspect  angle.  If  the 
value  of  o  is  near  0°  or  180°,  then  the  two  lines  are  nearly  parallel, 
and  the  distance  between  them  will  remain  unchanged  for  relatively 
large  changes  in  the  aspect  angle.  If  the  value  of  0  is  near  +90° 

(right  angles)  then,  the  two  lines  are  crossing  each  other,  and  relatively 
small  variations  in  aspect  angle  will  result  in  larger  values  of  D. 

Thus  radar  data  which  has  values  of  0  near  ±90°  will  have  better  sta¬ 
tistical  properties  for  target  identification  than  radar  data  with 
values  of  0  near  0°  or  180°. 

The  data  used  in  the  nearest  neighbor  example  (Missouri  (1:700 
and  Sverdlov)  was  also  used  to  compute  D  and  0  for  groups  of  harmonics 
over  limited  aspect  angle  ranges  as  shown  in  Table  4.  In  this  table, 
groups  of  4  harmonics  were  used  over  aspect  angle  scans  of  the  broad¬ 
side,  how  and  stern.  As  can  be  seen,  this  table  indicates  that  the 
best,  frequency  range  (largest  value  of  0)  occurs  for  harmonics  4  through 
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Dist.  (min.)  and  Skew  Angle  in  Hyperspace. 
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7  at  broadside,  harmonics  6  through  9  at  bow  on,  and  for  harmonics 
?  through  5  or  5  through  8  (a  tie)  for  the  stern  view.  These  harmonics 
are  preferred  with  more  than  a  two  to  one  margin  for  the  bow  or  stern 
aspects.  The  broadside  aspect  choice  of  frequency  range  is  less  clear 
if  only  the  value  of  D  is  considered.  If  the  value  of  theta  is  included 
in  the  argument,  however,  the  choice  of  the  4  to  7  harmonic  range  for 
optimum  target  identification  is  strongly  indicated.  The  choice  of 
frequency  range  for  the  stern  aspect  is  complicated  by  the  introduction 
of  o.  In  this  case,  large  values  of  D  are  accompanied  by  values  of  8 
which  are  less  desirable  ( i . e . ,  close  to  0°). 

An  important  use  of  this  computation  would  be  in  the  compari¬ 
son  of  the  value  of  0  to  the  expected  measurement  errors.  If  D  is 
significantly  larger  than  the  expected  measurement  errors,  then  the 
data  from  the  OTH  radar  return  from  ships  as  studied  here  would  be  in¬ 
herently  useful  for  ship  identification.  Since  the  error  statistics 
of  radar  cross  section  measurements  by  full  scale  systems  is  not  yet 
available,  this  comparison  is  not  yet  possible.  A  comparison  with 
the  average  magnitude  of  the  backscatter  return  is  possible,  however. 

The  results  are  given  below  for  the  values  used  in  Table  4. 

broadside  D  =  16.0,  Amp  (typical)  =  40. 
ratio  =  0.4 

stern  0  =  28.2  AMP  (typical  =  10 
ratio  =  2.82 

bow  D  =  62.1,  Amp  (typical)  =  10 
ratio  =  6.21 

These  results  are  very  encouraging.  The  implication  is  that  errors 
on  the  order  of  3  (stern)  or  6  (bow)  times  the  typical  backscatter 
measurement  level  can  be  tolerated  and  the  data  sets  will  still  be 
separated  in  terms  of  the  features  available  for  identification.  The 


broadside  results  indicate  that  the  cross  section  must  be  measured 
within  40%  of  the  value  of  a  typical  magnitude  in  orde>~  to  insure 
feature  set  separation.  The  magnitudes,  on  the  other  hand,  will  oe 
larger  and  the  absolute  measurement  accuracy  required  at  broadside 
may  be  the  same  as  at  bow  or  stern  aspects.  These  results  will  be 
discussed  further  in  the  last  section  of  this  report. 

V.  TIME  DOMAIN  FEATURES 

This  section  discusses  the  application  of  Transient  Electromagnetic 
Signature  research  to  the  problem  of  identifying  distant  ships  and 
aircraft  from  Over-the-Hori zon  (OTH)  Radar  returns.  First,  the  signature 
quantities  are  defined,  and  some  past  research  on  the  signatures,  th. 
nature  and  their  utilization  is  discussed.  Next  recommended  reouire"ients 
for  data  collection  based  on  the  nature  of  target  transient  signataa-es 
are  discussed.  Then  two  important  approaches  to  ship  identification 
are  described  assuming  that  the  necessary  OTH  data  are  available. 

The  first  approach  involves  the  information  contained  in  the  initial 
portion  of  the  transient  echo  return  waveform,  and  infers  geometrical 
quantities  of  the  body.  The  second  approach  involves  the  latter  portion 
■if  the  time-domain  signature  and  relates  to  the  characteristic  resonances 
of  fhe  object.  A  final  section  briefly  reviews  the  application  of 
those  approaches  to  over-the-hori zon  aircraft  identification. 

The  use  of  the  transient  radar  signature  of  object  was  first 
discussed  by  Kennaugh  and  Cosgriff  in  1957  [51.  Since  then,  considerable 
research  effort  has  been  done  concerning  basic  transient  signature 
properties  [61,  calculation  of  transient  signatures  for  specific  bodies 
[71,  and  antennas  for  possible  transient  radar  use  [8].  The  time-  lorn. a  in 
transient  signature  used  in  this  effort  is  the  ramp  response  wavdOr  , 
which  was  first  suggested  for  radar  identification  hy  Kennaugh  and 
Mol  I  at t  (ol.  This  signature  is  defined  as  the  far-zone  backseat ! or. >d 
lime  domain  waveform  of  a  target  i  1  luminal."  I  by  a  traveling  plninr 
transverse  electromagnetic  field  unit  slope  discontinuity.  The  impul,1, 
top,  and  ramp  response  signatures  of  a  sphere  are  shown  in  Figure  7. 


The  ramp  response  signature  is  the  second  integral  of  the  impulse, 
response  ol  a  target,  and  hence  shares  several  of  its  useful  properti.-s. 

I  he  ramp  response  is  unique  with  respect  to  target  shape,  orientation, 
and  material  composition.  It  is  the  inverse  Fourier  transform  of  a 
target's  complex  backseattering  frequency  spectrum  multiplied  by  the 
factor  (l/(j<o))2: 

R( t)  (i.)  (i) 

2 

Because  of  the  (l/(jw))  factor,  Kennaugh  and  Moffat*  predicted 
that  satisfactory  approximate  ramp  response  signatures  could  be  obtained 
with  narrower  bandwidth  interrogating  signals  than  for  the  impulse 
response.  For  example,  ten-harmonic  approximations  of  the  sphere  impulse 
and  ramp  response  waveforms  are  shown  in  Figure  4.  The  fundamental 
frequency  has  wavelength  approximately  five  times  the  sphere  diameter. 

As  can  be  seen,  the  major  features  of  the  ramp  response  waveform  are 
faithfully  reproduced. 

It  should  be  noted  that  these  signatures  are  all  functions  of 
aspect  angle  and  polarization  of  the  interrogating  signal  with  reseed 
to  the  target.  Thus,  quanti tati vely 

gipA^iM  =  f  R(p,e,<j.,t)e"J’ut  dt 

(jo))2  o 


whore 


2 

TTC  | G (  p,  0,<f>,  jw)  |  5  o(p,0,<j>,ju>) 

and  p,0,<t  emphasize  the  dependence  on  polarization  and  interrogation 
iqriul  direct  ion. 
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Figure  4.  Ten-harmonic  approximate  sphere  transient  signatures. 


20 


Data  Regu ir (.merits  for  Generating 
Ramp  Response  Signatures 

Because  of  the  Fourier  Transform  relationship  it  is  obvious  that 
signature  data  may  be  gathered  either  by  the  use  of  time-domain  transien 
interrogating  radar  signals,  or  by  measuring  the  complex  backscattered 
cross-section  of  a  target  at  sampled  frequencies  in  the  frequency  domain 
As  was  already  shown,  adequate  ramp  response  signatures  for  simple 
shapes  are  obtained  from  data  at  as  few  as  10  harmonically  related 
frequencies. 


Based  on  known  ramp  response  characteri sties,  we  can  calculate 
the  frequency  sampling  interval  for  valid  signature  reconstruction, 
and  the  relative  importance  of  the  different  frequencies.  First,  it 
is  known  that  the  ramp  response  signature  is  of  finite  duration.  Using 
a  conservative  estimate  based  on  past  measurements,  the  duration  of 
the  signature  is  less  than  6  transit  times  across  the  length  of  the 
object,  i.e., 


R(p,M,t-t0) 


0  (t-tQ)  <  0 


0  (t-tQ) 


> 


c 


where  this  ramp  response  in  retarded  time  begins  at  time  tQ,  &  equals 
object  length,  and  c  is  the  speed  of  light. 

For  a  time-limited  function,  application  of  the  Shannon  samplino 
theorem  produces  a  required  frequency  sampling  interval.  Specifically, 
given 


f(t)  =  0  for  | t |  >  T, 


then 


F(<n) 


sin(g)T-mr) 

(jT-nir 


?1 


Thus,  for  a  given  object  length  l,  the  scattering  signature  of  the 
object  should  be  sampled  at  intervals  of 


c  nn 
n  ~  TT2?) 


and  T  =  \  (f^)  ; 


o 

nc  n  x  3x10  meters/sec 

vn  "  6l  "  6  S>  , 

meters 


This  formula  implies  that  for  an  aircraft  carrier  of  300  meters  length, 
the  sampling  frequency  interval  should  be  183  KHz.  For  a  destroyer 
of  125  meters  length,  a  sampling  frequency  interval  of  400  KHz  or  less 
is  sufficient.  It  should  be  stressed  that  these  numbers  are  derived 
from  worst-case  signature  characteristics  and  are  sufficient  for  all 
cases.  For  some  aspects  and  interrogating  polarizations,  sampling 
this  often  in  frequency  may  not  be  necessary.  On  the  other  hand,  main¬ 
taining  such  sampling  intervals  in  these  cases  provides  some  redundant 
information,  in  effect  improving  the  signal  to  clutter  performance 
of  the  signature  measurement  system. 


Another  important  feature  of  the  ramp  response  signature  is  its 
| l/m‘ |  weighting  of  frequency  domain  cross-section  data.  This  weighting 
emphasizes  the  low-frequency  scattering  characteristics  of  targets. 

For  example,  it  means  that  acceptable  ramp  response  signatures  may 
be  derived  from  frequencies  below  12  MHz  for  all  ships  larger  than 
125  meters  long. 


Target  Geometry  vs.  Ramp  Response  Signature 

One  of  the  most  attractive  features  of  using  ramp  response  sig¬ 
natures  derived  from  over-the-horizon  radar  is  that  the  first  portion 
of  the  ramp  response  approximates  the  target  "profile  function"  for 
metallic  targets  [9].  A  calibrated  ramp  response,  containing  both 
ampl  itude  and  phase  vs.  frequency  input,  is  approximately  equal  to 
the  target  cross-sectional  area  vs.  distance  along  the  line  of  sight. 
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Thus,  particularly  for  near-bow  and  near-stern  look  angles,  it  is  seen 
that  such  a  signature  carries  considerable  information  on  the  gross 
features  of  a  ship. 

The  above  relationship  has  been  demonstrated  using  measured  data 
for  simple  shapes  like  cones,  cylinders,  spheroids  etc.  [9, 10").  It 
was  examined  for  ship  targets  on  this  program,  using  the  measured  model 
scattering  data  (Table  1). 

The  size  of  all  the  measured  ship  models  were  too  large  for  the 
frequency  sampling  criterion  just  discussed  to  be  satisfied.  Therefore, 
aliasing  is  present  in  all  of  the  reconstructed  ramp  response  waveforms. 
The  smallest  ship  models  were  chosen  for  study,  because  it  was  expected 
that  their  data  would  be  most  accurate.  The  amplitude  spectra  for 
these  two  ships  for  nose-on  incidence  are  shown  in  Figure  5  (Hayanami) 
and  Figure  6  (Shimokaze)  [3 1 .  It  is  very  interesting  to  note  that, 
despite  the  fact  that  the  two  destroyer  models  are  almost  of  identical 
size  and  have  many  similarities  in  their  gross  structural  outlines, 
there  is  a  relative  null  feature  in  the  Hayanami  spectrum  which  does 
not  appear  in  the  Shimokaze  spectrum.  There  is  not  a  data  point  to 
define  the  actual  depth  of  this  null  which  was  generated  by  the  spectrum 
smoothing  process.  However,  some  features  of  the  ramp  responses  tend 
to  confirm  its  existence  and  explain  why  it  happens. 

The  approximate  ramp  response  waveforms  for  these  ships  are  shown 
in  Figure  7,  along  with  correctly  scaled  profile  outlines  of  the  two 
ships.  It  is  noted  that  there  appears  to  be  experimental  model  position¬ 
ing  error  in  the  nose-on  look  angle  measurement  of  the  Shimokaze. 

It  should  probably  be  moved  -  3  cm  forward  to  agree  with  the  stern 
look  angle  waveform. 

There  is  an  important  difference  in  the  waveforms  of  these  two 
'.flips,  which  probably  explains  the  relative  null.  It  is  seen  that 
the  Hayanami  has  two  peaks,  corresponding  to  the  forward  and  aft  vertical 


Magnitude  spectrum  vs.  scaled  frequency  for  the 
Hayanami  (1:700)  -  bow  on. 


s 

SCALED  FREQ.  (H«) 


SHIMOKAZE  (1/700) -BOW  ON 

Magnitude  spectrum  vs.  scaled  frequency  for  the 
Shimokaze  (1:700)  -  bow  on. 
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Figure  7.  Ramp  response  results  and  ship  outlines. 
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masts.  The  relative  null  in  the  Hayanami  frequency  sepctrum  simply 
reflects  an  interference  effect  between  the  approximately  equal  strength 
echos  of  the  fore  and  aft  masts.  The  Shimokaze  does  not  have  an  aft 
mast  structure,  and  hence  has  only  one  major  peak  in  its  ramp  response. 

Two  conclusions  are  drawn  from  this  information.  First,  the 
measured  ramp  response  data  shows  promise  as  a  tool  for  determining 
which  structures  on  the  ship  make  major  contributions  to  the  scattering 
signature  vs.  frequency,  time-delay,  and  aspect  angle.  Thus,  even 
though  it  is  probably  impractical  to  reconstruct  a  ramp  response  signature 
with  an  operational  system,  it  would  be  important  to  perform  further 
rainp  response  research  on  models  in  order  to  guide  a  selection  of  re¬ 
liable,  geometrical ly  meaningful  features  in  the  frequency  domain  for 
the  operational  system.  A  second  conclusion  is  that  for  these  two 
ships,  the  scattering  at  near  on-axis  look  angles  primarily  reflects 
the  vertical  masts. 

VI.  SHIP  TARGET  CHARACTERISTIC  RESONANCES 

Concerning  the  complete  transient  response  (impulse,  step,  or 
rainp)  of  a  ship,  the  time  domain  waveform  contains  the  initial  profile 
function  just  discussed  followed  by  a  decaying  complex  exponential 
portion  which  is  characterized  by  the  natural  resonances  of  the  object. 

It  is  important  to  note  that  while  the  amplitudes  and  phase  of  their 
excitations  may  change  as  a  function  of  direction  and  polarization, 
the  poles  in  the  complex  plane  which  correspond  to  the  natural  resonances 
are  invariant. 

Two  theoretical  and  experimental  efforts  using  natural  resonance 
for  target  identification  have  been  studied  and  evaluated  for  potent,  i  a  1 
application  to  OTH  radar  returns.  The  first  effort  involved  transimit 
radar  measurements  of  buried  metallic  and  non-metallic  targets  flP. 

The  second  effort  is  the  one  whose  model  measurement  data  has  boon 
used  in  this  report,  concentrating  on  radar  frequencies  of  1.4?,  ?.M4 
...  14.?  MHz  incident  on  ships. 
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The  underground  radar  study  dealt  with  metallic  and  non-metallic 
objects,  few  of  which  were  strongly  resonant.  A  target  set  is  shown 
in  Figure  8,  and  a  set  of  poles  obtained  from  measured  data  is  shown 
in  Figure  9  for  one  target.  These  data  show  that  for  the  first  few 
poles,  their  location  is  almost  excitation  invariant.  It  is  also  im¬ 
portant  that  for  these  targets  which  are  not  strongly  resonant,  the 
poles  can  be  extracted  from  measured  scattering  data. 

Poles  have  also  been  extracted  for  the  ship  models  [3],  although 
their  locations  seem  to  be  somewhat  dependent  on  ship  orientation.  There 
are  two  possible  reasons  for  this:  1)  The  data  have  aliasing,  and 
2)  the  ship  models  are  more  complex,  and  therefore  should  have  more 
natural  resonances  than  the  simple  underground  targets. 

Techniques  have  been  devised  for  target  identification  using 
target  poles.  It  is  found  that  identification  requires  only  a  portion 
of  the  total  excitation  spectrum,  and  is  independent  of  target  aspect 
angle.  Based  upon  these  features,  the  characteristic  resonances  of 
ships  are  regarded  as  promising  candidates  for  use  as  features  for 
identification  using  OTH  data.  These  resonances  need  to  be  studied 
further,  with  a  data  base  that  is  not  degraded  by  aliasing. 

Finally,  previous  studies  have  shown  that  transient  signatures 
are  a  useful  tool  for  guiding  the  selection  of  features  for  aircraft 
identification.  In  particular,  identification  techniques  using  fre¬ 
quencies  pertinent  to  OTH  radar  have  already  been  developed  [12]. 

VII.  CONCLUSIONS  AND  OVERVIEW 

A  number  of  studies  have  been  performed  with  the  goal  of  develop¬ 
ing  a  system  for  the  automatic  identification  of  ships  using  the  radar 
return  from  an  H.F.  over  the  horizon  (OTH)  radar  system.  This  phase 
of  the  research  has  been  involved  with  using  available  techniques  and 
an  available  data  base  to  study  the  general  identifiability  of  ships 
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Figure  8.  Physical  characteristics  of  the  subsurface  targets 
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mine- 1  ike  target  at  different  antenna 
locations  in  icy  ground. 


based  on  OTH  radar  returns.  The  results  are  quite  encouraging.  It 
has  been  shown  that  in  general,  ships  are  inherently  identifiable  from 
il.F.  OTH  radar  data.  The  features  cluster  and  separate  far  enough 
to  accommodate  large  error  margins.  This  means  that  the  application 
of  specific  pattern  recognition  techniques  to  this  problem  can  be  expected 
to  yield  good  results.  The  choice  of  the  particular  algorithms  to 
employ  for  automatic  target  identification  will  depend  on  further 
studies  of  data  and  of  the  error  mechanisms  and  statistics  found  in 
operational  OTH  systems. 

Studies  in  the  time  domain  complement  these  results.  These 
studies  have  shown  that  the  large  significant  structures  on  a  ship 
are  the  major  contributors  to  the  H.F.  radar  return.  Quantitative 
values  of  the  sampling  density  in  the  frequency  domain  were  derived 
using  time  domain  techniques. 

The  frequency  domain  and  the  time  domain  results  are  also  in 
agreement  with  regard  to  the  aspect  angles  which  are  important  for 
future  study.  Both  techniques  indicate  that  broadside  data  from  a 
ship  is  most  difficult  to  use  for  target  identification.  One  reason 
suggested  for  this  is  that  the  broadside  data  are  stronger  than  the 
bow  data  because  of  the  large  return  from  the  hull  of  the  ship.  Since 
the  identification  characteristics  have  been  shown  to  be  related  to 
the  arrangement  of  the  superstructure,  the  return  from  the  hul 1  of 
the  ship  at  broadside  can  be  thought  of  as  a  clutter  contribution. 

This  "clutter  term"  is  large  at  broadside,  but  small  at  bow  and  stern 
aspects.  The  return  from  the  superstructure,  on  the  other  hand,  is 
of  the  same  order  of  magnitude  for  the  bow  and  broadside  cases.  Studies 
in  this  area  should  be  continued. 

Future  studies  of  H.F.  OTH  radar  identification  of  ships  hav» 
a  prerequisite  for  a  specialized  data  base.  The  data  base  should 
be  designed  for  ship  target  identification  using  data  from  H.F.  radar. 

This  means  that  a  set  of  data  for  full  aspect  angle  coverage  and  0° 
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to  30°  elevation  angle  coverage  for  a  number  of  ships  should  be  obtained. 
The  frequency  coverage  should  be  consistent  with  the  constraints  mentioned 
in  Section  V  of  this  report.  The  ships  to  be  studied  should  be  chosen 
based  on  the  need  to  support  a  full  scale  OTH  H.F.  radar  system.  This 
data  base  should  be  of  sufficient  quantity  to  allow  development  and 
testing  of  the  performance  of  target  identification  algorithms  relevant 
to  the  full  scale  operational  system.  Comparisons  between  the  model 
measurements  and  the  full  scale  measurements  will  be  made.  Testing 
of  target  identification  algorithms  on  the  available  full  scale  measure¬ 
ments  can  then  begin. 
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APPENDIX  A 


The  raw  data  used  in  the  target  identification  studies  presented 
in  this  report  is  given  here.  The  raw  data  was  in  the  form  of  amplitude 
and  phase  measurements  of  the  signal  backscattered  from  a  scale  model 
ship.  The  ship  models  used  were  the  Sverdlov  and  the  Missouri  (1:700 
scale).  The  amplitude  (measured  in  cm)  is  given  as  the  square  root 
of  the  radar  cross  section.  The  phase  is  given  as  degrees  absolute 
with  a  phase  reference  determined  by  a  reference  sphere.  The  reference 
point  for  the  phase  data  is  a  point  at  the  center  of  the  hull  of  the 
ship  at  the  waterline.  The  Harrington  phase  reference  convention 
(increasing  distance  along  the  ray  path  implys  decreasing  phase  values) 
is  used  here. 

The  values  of  amplitude  and  phase  are  presented  as  polar  plots 
so  that  comparisons  between  data  from  various  ships  or  various  aspect 
angles  may  be  made.  The  zero  phase  direction  is  along  the  x-axis  (to 
the  right) . 

The  data  are  given  in  three  groups;  bow  data,  stern  data,  and 
broadside  data.  The  polar  plots  are  identified  by  a  code  of  th*>  form 
MISxxx  or  SVExxx.  The  MIS  and  the  SVE  part  of  the  identification  stand 
for  Missouri  and  Sverdlov  respectively.  The  xxx  gives  the  aspect  angle 
measured  in  degrees  from  bow  on. 
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APPENDIX  B 


The  following  tabulations  are  the  nearest  neighbor  computation 
results  from  this  report.  They  were  summarized  earlier  in  this  report. 

The  columns  are  headed  with  the  code  MISxxx  or  SVExxx.  This  represents 
either  the  Missouri  or  the  Sverdlov  at  aspect  angle  xxx  (given  in  degrees) 
The  distance  (nearest  neighbor)  from  the  data  for  the  ship  given  in 
the  first  column  to  the  data  for  the  ship  given  at  the  top  of  the  fol¬ 
lowing  columns  is  printed  at  the  points  of  intersection.  As  can  be 
seen,  the  distance  for  identical  ship  data  values  is  always  zero. 

The  remaining  numbers  are  the  result  of  the  nearest  neighbor  distance 
computation.  No  normal ization  has  been  done.  The  minimum  value  in 
the  horizontal  row  is  underlined.  It  is  underlined  once  if  the  minimum 
is  for  the  same  ship  and  twice  if  the  minimim  is  for  the  other  ship. 

In  the  upper  left  are  two  numbers.  The  first  number  is  the  number 
of  the  first  harmonic  used  in  the  analysis,  and  the  second  number  is 
the  number  of  the  last  harmonic  used.  The  numbers  are  followed  by 
either  a  T  or  an  F.  The  T  indicates  that  phase  data  has  been  included 
in  the  computations,  and  the  F  indicates  that  it  has  not. 

The  algorithms  used  for  these  tabulations  is  given  below. 

For  amplitude  data  only: 
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j-p  J  J 


where: 

x.  and  y.  are  amplitudes  of  the  jth  harmonic  of  two 
J  J  data  sets  (x  and  y) 

p  is  the  starting  harmonic  number 

q  is  the  final  harmonic  number 

2  2 
D  is  tabulated  since  a  minimum  of  D  implys  a  minimum 

of  D. 


for  amplitude  and  phase  data: 


q  2  q  2 

D2  =  Z  Ursine  -y.sine  )  +  £  (x.cose  -y.cose  ) 

j=P  J  xj  J  y3  j=P  J  J  3  yJ 

where: 

o  and  o  are  the  phase  angles  of  the  jth  harmonic  of 
xj  yj  two  data  sets  (x  and  y). 
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APPENDIX  C 

HYPERLINE  DERIVATION 


Data  from  a  frequency  scan  on  a  target  defines  a  vector*  in 
multidimensional  space. 

V’(°0  «  xpl(ot)^1  +  xp,2(ot)i2  +  +  xpn(ot)in 

pth  target 

n  components  (or  dimensions) 

(1  or  2  components  per  frequency  in  the  case  of 
the  ship  data) 

V^(a)  will  trace  a  line.  Assume  a  fit  to  a  straight  line  is  meaningful. 

y«>  -  v  *  bp 

where: 

mp  -  "pi  h  *  V  '2  *  •••  +  V 

^  A  A 

=  b  .  ii  +  i9  +  ...  +  i 
p  pi  I  p2  2  pn  n 

The  distance  between  two  points  on  two  lines: 

dl,2^“l’a2^  =  I vi 
The  skew  angle  <|>j  2: 

M,-M2 

cos  *1,2  ’  IW,1-|H2| 

If  Vj  and  V2  are  "good"  estimates  of  V^  and  V2,  then  d^  and  4>-j 2 
are  measures  of  the  classification  potential  of  versus  V£. 


*A1 1  vectors  in  this  appendix  are  capatalized.  All  scalars  are  small 
letters. 
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| D |  =  d  is  minimum  when: 

D-M1  =  0  and  D-Mp  =  0 
D  -  VrV2 
V  =  Mch-B 

(VrV2)-M1  =  V1*M1  -  V2-M1  =  0 
Thus:  [Scalars] 

(M-j  *M-| ) otj  -  (M2-M])o(2  +  (BrB2)  *M1  *  0 

(M]-M2)o1  -  (M2-M2)o2  +  (BrB2)*M2  =  0 

At  this  point,  there  are  two  scalar  equations  in  two  unknowns  ( 
SO:  Find  aj  and  «p 

Gives  d  minimum  and  4>. 


^  and  m2). 
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